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ABSTRACT 

As the most common environment in the universe, groups of galaxies are likely to contain a significant 
fraction of the missing baryons in the form of intergalactic gas. The density of this gas is an important 
factor in whether ram pressure stripping and strangulation affect the evolution of galaxies in these 
systems. We present a method for measuring the density of intergalactic gas using bent-double radio 
sources that is independent of temperature, making it complementary to current absorption line 
measurements. We use this method to probe intergalactic gas in two different environments: inside a 
small group of galaxies as well as outside of a larger group at a 2 Mpc radius and measure total gas 
densities of 4±1^ x 10™ 3 cm -3 and 9±3j^° x 10~ 4 cm -3 (random and systematic errors) respectively. 
We use X-ray data to place an upper limit of 2 x 10 6 K on the temperature of the intragroup gas in 
the small group. 

Subject headings: intergalactic medium - galaxies; jets - galaxies: clusters: general 



6 



m 

o 

oo 

O 

X 



1. INTRODUCTION 

In the local universe, attempts to measure the baryon 
content in the form of stars, hot X-ray emitting gas, and 
cold gas account for only one third of the baryo n density 
seen at high redshift (|Fukugita fc Pee bles 2004) . Simula- 
tions predict that the remaining "missing baryons" exist 
in a warm-hot intergalactic medium (WHIM) that was 
shock heated during large-scale structur e formation to a 
temperature range of 1 5 < T < 10 7 K (|Cen k. Ostrikerl 
l2006t iDave et al.l l200lh . Observational confirmation of 
the WHIM comes primarily from UV absorption lines i n 
the spectra of low redshift quasars ( Tripp et al.l [2000) . 
However, UV absorption line observations are only able 
to probe a limited temperature range of WHIM gas and 
there is continued debate as to whether the absorption 
originates in the extended gaseous halos of individual 
galaxies, the int ergalactic medium (IGM), or from large- 
scale filaments (| Savage et al.l [20031 ). 

Galaxy groups, the most common environment in the 
local universe, likely contain a significant fraction of 
the total baryonic ma ss within their intragroup medium 
(jFukugita et alj|1998f ). Until now the only probes of this 
gas have been UV absorption line and X-ray observa- 
tions. In the few groups that have a hot intragroup 
medium X-ray observations can determine the emission 
measure, but without some knowledge of the geometry 
the density of this gas cannot be disentangled. The X-ray 
surface brightness typically falls below detectable levels 
at radii much smaller than the virial radius, while the 
implied gas mass i s still increasing at least linearly with 
radius (jMulchaevi [2000) . Groups with a bright quasar 
behind them have been observed in an attempt to detect 
absorption from the intragroup medium. In cases where 
intervening gas is detected, assumptions about the ex- 
tent of the absorbing system, the metallicity, ionization 
fraction, and the intensity of the ionizing radiation field 
lead to estimates of the intragroup medium density on 
the order of n ~ 10 -4 — 10 
ITripp et al.lll99"i . 



5 cm" 3 (jPisano et all 120041 : 



Additionally, galaxy groups, according to the hier- 
archical scenario of the formation of large-scale struc- 
ture, are the building blocks of rich clusters of galax- 
ies. Groups are important sites in which to investigate 
the physical mechanisms responsible for the observed 
morphology and SF R-density relations ([Goto et alj|2003t 
iGomez et al.ll2003h . However, the effectiveness of mech- 
anisms like ram pressure stripping and strangulation on 
galaxy evolution in the group environment are not well 
understood due to the uncertainty in the density of the 
intragroup medium. Some observational evidence does 
exist for the stripping of hot gas from galaxies by the 
intragrou p medium, mainly in features like X-ray sho cks 
and tails (ISivakoff e t al. 2004; Rasmu ssen et al.ll2006l ). as 
well as HI observations that show com pression features 
opposite extensions in cold d isk gas (|Kantharia et al.l 
[20051 iMcConnachie et ai]l20"ol . 

Another way to measure the density of the IGM is 
through its effect on the jets of radio galaxies. As a dou- 
ble lobed radio source travels through intergalactic gas 
ram pressur e caus es its jets to be swept back. Prior to 
iBurns et alJ (|1987[ ) it was believed that the combination 
of the dense intracluster medium and large galaxy veloci- 
ties needed to form these bent-double radio sources could 
only be found in rich clusters. While many bent-doubles 
do reside in dense clusters, a surprising number are 



fou nd in lower mass environments like groups of galax- 
[19951: iBlanton et al.l[200l . 



ies (|I^rs_et~al . 1978; Venk atesan et al.iri994t iDoe et alJ 
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We present measurements of the density of intergalac- 
tic gas from radio, optical, and X-ray observations of two 
bent-double radio sources and their sur rounding envi- 
ronme nts. Both sources were identified by IBlanton et al.l 
(|2001l) as existing in environments with richness less than 
Abell class 0. We adopt H a = 73 km s^ 1 Mpc" 1 . 

2. DATA 

2.1. Radio Observations 

We observed the FRI radio galaxies FIRST J124942.2+ 
303838 and SDSS J085331.86 + 232400.0 (hereafter SI 
ev ^§^?S^m i$Sy) in March of 2006 with the Giant Me- 
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trewave Radio Telescope (GMRT). The GMRT consists 
of thirty 45 meter diameter antennas with half arranged 
in a compact 1 km diameter, randomly distributed, cen- 
tral square and the remaining in three extended arms 
with fixed baselines ranging from 0.1 to 25 km in length. 
The observations were dual frequency 610/235 MHz us- 
ing both the upper and lower sidebands. In this paper 
we will discuss only the single polarization 610 MHz data 
which have an effective bandwidth of 12 MHz in each 
sideband. We observed a flux calibrator (3C286 or 3C48) 
at the beginning and end of each observation and every 
40 minutes we observed a phase calibrator for 6 min- 
utes. We have a total of 5.5 hours integration on each 
source. The data were calibrated in AIPS using standard 
tasks. Imaging and self-calibration were done with mul- 
tiple facets across the 0.7° (half-power beam width) field 
of view. 

2.2. Optical Observations 

Multi-object spectroscopy was performed in May 2007 
and April 2008 using HYDRA on the WIYN 3.5m tele- 
scope on Kitt Peak for the source SI. Fibers from the 
blue cables were placed on galaxies chosen from Sloan 
Digital Sky Survey (SDSS) photometric redshifts to have 
a redshift similar to the radio source, and additional 
fibers were placed on blank sky positions. The 600@10 
Zepf grating was used in first order to give spectra with a 
dispersion of 1.4 A per pixel covering a wavelength range 
from 4600 to 7200 A. CuAr lamp calibration spectra 
were used for the wavelength calibration and an averaged 
sky spectrum was subtracted from each object spectrum. 
Spectra were obtained for ~ 60 galaxies overall. The 
data were reduced in IRAF using the dohydra package. 
Cross correlation of object spectra with g alaxy template 
spect ra was performed using the RVSAO ([Kurtz fc Mink! 
1 1998( 1 package in IRAF. 

2.3. X-ray Observations 

We obtained archival Chandra observations of both 
sources taken originally in January and February 2005 
(PI Blanton). Exposure times were 35.17 and 47.19 ksec 
for SI and S2, respectively, and both sources were cen- 
tered on the back illuminated ACIS-S3 chip. 

The event files were filtered to include events with en- 
ergies between 0.3 keV and 7.0 keV. Point sources were 
identified using task wavdetect and removed. Images 
were made by smoothing the data in two dimensions 
with a gaussian extending to 5<r and with a full-width 
at half maximum of 6 pixels on both axes. Images were 
also made using an adaptive smoothing with a minimal 
significance signal to noise ratio of 3. Spectra were ex- 
tracted from a region 300 kpc in radius around the radio 
source as well as a background region. These sources 
are distant enough that they fit on the back illuminated 
ACIS-S3 chip with area left for a generous background 
region. 

3. METHOD 

For a relativistic, hydrodynamic jet, the time- 
independent Euler equation describes t he balance of in- 
terna l and external pressure gradients (iBegelman et al.l 
H97<l Uones fe OweiJll979t iBurns fe Owenlll980D . 
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where Pi GM v 2 , is the external ram pressure felt by the 

radio galaxy as it travels through the IGM, wT 2 /3 2 is 
the relativistic enthalpy density inside the jet, h is the 
width of the jet, and R is the radius of curvature of the 
jet. In our data the jet widths are not resolved and are 
set by the beam size of the individual observations. The 
radius of curvature, R, is found by fitting a circle by 
eye along the jets through the core. This radius is il- 
lustrated in the left panels of Figure [T] The velocity 
of the radio galaxy, v ga i , is estimated using the velocity 

dispersion of the group as y3a v . The enthalpy density 
is written as w = e + p, where e is the energy den- 
sity and p is the internal pressure. For the jets whose 
particle population is ultrarelativistic, e = 3p + pj e tc 2 
and thus w = AP m in + Pjetc 2 assuming that the internal 
pressure in the jets is dominated by the minimum syn- 
chrotron pressure, P m in- While Pjetc 2 is not well known, 
it is anticipated to be much less than 4P m ,;„ and we d o 
not include it in our calculations ([Bridle fc Per levll 19841) . 
We adopt (3 = v/c — 0.54 ± 0.18 as the distribution of 
obse rved jet speeds in FRI rad io galaxies with straight 
jets (jArshakian fc LongaiHl2004l ) which is co nsistent with 
speeds in wide-angle tailed radio galaxies (|Jetha et al.l 
120061 ). The relativistic factor T is the usual (1 - /3 2 )" 1/2 . 

The internal pressure is calculated using the mini- 
mum synchrotron pressure as outlined in lO'Dea fc Owenl 
(fl987h . 

Pmin = (2rr)-* (j^J [ci 2 L rad (l+fc)( ( /)F)- 1 ]? ergs cm" 3 

(2) 

where c\i is a constant th at depends on the spectral in- 
dex and frequency cutoffs (|Pacholczvklll970T ). k is the ra- 
tio of relativistic proton to relativistic electron energy, <f> 
is the volume filling factor, V the source volume, and L ra d 
the radio luminosity of the jet. We measure this internal 
pressure at the position immediately before the jet bends, 
always excluding the core. Standard equipartition of en- 
ergy between relativistic particles and magnetic fields is 
assumed. We follow the general assumptions that fc = 1, 
(f> = 1, and the jets are cylinders uniformly filled with 
magnetic field s and relativistic p articles. We use VLA 
FIRST survey (|Becker et al.lll995[ ) data at 1420 MHz and 
our GMRT 610 MHz data to determine the spectral in- 
dex of the synchrotron emission. Our values for the min- 
imum synchrotron pressure in the jets are in good agree- 
ment with other measurements for similar radio sou rces 
([Venkatesan et al.lll994 IWorrall fc Birkinshawil2000f ). 

4. GROUP REDSHIFTS AND VELOCITY 
DISPERSIONS 

In order to estimate the velocity of the radio galaxy 
through the IGM we use the optical data to examine the 
environments of these sources. Small group memberships 
and instrumentation limit our ability to obtain large ra- 
dial velocity samples. Additionally, structures associ- 
ated with these radio sources may not have achieved the 
dynamical equilibrium that would result in a Gaussian 
distribution of radial velocities among members. With 
these considerations in mind we use the robust statistical 
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biweight estimators for location and scale (|Beers et al.l 
1990) to determine the group redshift and velocity dis- 
persion for each sample. These estimators do not assume 
a Gaussian parent population and are not strongly influ- 
enced by outliers which is especially important for small 
samples. 

We combined redshift information from our WIYN 
spectroscopy (only for SI), as well as spectroscopic and 
photo metric redshifts from the Sloan Digital Sky Survey 
DR6 (jAdelman-McCarthv et~al"1l2008l ) . We use only pho- 
tometric redshifts whose errors are less than 5% of their 
value for SI and 10% for S2. 

We calculate a peculiar velocity for a galaxy with red- 
shift z in the rest frame of the group using 

Vpec C\Z Zgroup)/(\ ~t~ Zgroup) • (3) 

The velocity dispersion, a group , will be taken as the dis- 
persion of the v pec values. We start by examining the 
redshift histogram for galaxies within a projected 6 Mpc 
distance from the radio source and the sky locations 
of galaxies within 2000 km s _1 . With so few redshifts 
available the final samples are chosen by eye based on 
the observed spatial clustering seen in Figure [T] which 
shows the spatial and redshift distributions of galaxies 
near the two radio sources. In the case of S2, the radio 
source appears offset spatially from the group. 10 of the 
galaxies circled in black have photometric redshifts and 
3 (including the radio source) are SDSS spectroscopic 
redshifts. We estimate the speed of this source through 
the IGM using its velocity relative to the mean redshift 
of the system it is approaching (z = 0.3081) and find a 
speed of 570 km s~ . The velocity of this source is un- 
certain although not unreasonable and we use a rough 
error of 10% of its value in future calculations. SI is 
associated with four other nearby galaxies, three have 
redshifts measured from our WIYN data and the fourth 
is a photometric redshift from SDSS. These five galaxies 
have a velocity dispersion of 250+^0 km s _1 . The 68% 
error bars on the velocity dispersion are determined by 
generating bootstrap samples using the final set of v pec 
values. 

5. IGM TEMPERATURE 

We use the X-ray data to examine t he temperature 
of the IGM in these systems. XSPEC (jArnaudl l200l 
is used to fit the spectra extracted from a region 300 
kpc in radius surrounding the radio source with point 
sources removed. We use Raymond models for thermal 
X-ray emission from a hot, optically thin gas, fixing the 
abundance at Z=0.1 and the normalization according to 
the measured number density shown in Table [TJ assum- 
ing it is constant throughout the volume. All models in- 
cluded a WABS component of Galactic absorption with 
values of 1.25 x 10 20 cm~ 2 and 3.23 x 10 20 cm~ 2 for 
SI and S2 respectively, b ased on the measurements of 
iDickev k Lockmanl (| 19901) . 

In the case of SI the number of counts in the region 
around the radio source after background subtraction 
is only 80 ± 35. The spectrum from this region was 
binned to give each measurement a 3er significance level 
above the background. A Raymond model was then com- 
pared with the data to determine the temperature the gas 
would need to produce counts at a 3er level above the cur- 
rent observations. Gas with a density of 4 x 10~ 3 cm~ 3 



and a temperature higher than 0.2 keV (2 x 10 6 K) should 
be observable with the current data. Thus, we put an up- 
per limit on the temperature of the IGM in this group 
at 2 x 10 6 K. 

After background subtraction we detect 94 ± 26 counts 
from the source region for S2. We do not expect the in- 
tergalactic gas around this source to be hot since it is 2 
Mpc from the center of a system of galaxies. The adap- 
tively smoothed image shows the X-ray counts concen- 
trated near the radio source. We interpret these X-rays 
as originating from Inverse-Compton scattering off the 
relativistic electrons in the radio source. The data are 
not adequate to distinguish between thermal and non- 
thermal models for the emission. Additionally, with the 
current data the density of this gas is too low to place an 
upper limit on the temperature were we to assume that 
these X-rays originated from a hot IGM. 

6. DISCUSSION 

We present two measurements of the density of inter- 
galactic gas using radio sources whose jets are bent back 
by ram pressure. SI is probing gas with a density of 
4± 1+2 x 10~ 27 g cm~ 3 (4± l+Jf X 10~ 3 cm -3 assuming 
a mean molecular weight of 0.6 and where the first set 
of errors are random and the second systematic) that is 
near the center of a small group of galaxies. In the case of 
S2 this source is probing gas that is a projected distance 
of 2 Mpc from the center of a system of galaxies with a 
density of 9±3t5°xl0- 28 g cm- 3 (9±3±5 5 xl0- 4 cm" 3 ). 
We have not accounted for thermal pressure in the jets, 
as well as our inability to resolve their true width. Ad- 
ditional thermal pressure or smaller jets will necessitate 
a larger IGM density to produce the observed jet geom- 
etry. If we have overestimated the radio galaxy speed 
or projection effects have led us to a smaller radius of 
curvature then the necessary density will decrease. 

The systematic errors reflect the range of IGM densi- 
ties possible given the range of jet speeds discussed in 
Section 3. These speeds are measured by assuming that 
flux asymmetries between the jet and counter-jet are the 
result of relativistic beaming effects. This beaming de- 
pends on the angle of the jets from our line of sight as well 
as the jet speed. By looking at large samples of objects 
and assuming a uniform distribution of angles, a range 
of jet speeds can be determined. We lack the inclination 
information nec essary to determine je t spee ds for indi- 
vidual sources. lArshakian k Longairl (|2004l ) point out 
that their distribution of jet speeds is asymmetric with 
the peak toward lower speeds and that selection effects 
may be pushing the mean toward higher speeds. Slower 
jet speeds would cause our IGM densities to tip towards 
the lower end of the range in values. 

The inclination of the source also affects the measured 
radius of curvature. We use a value which is uncorrected 
for inclination to calculate the IGM densities. SI appears 
relatively unprojected whereas S2 is less easily character- 
ized. If we assume that a Hubble flow distance is accurate 
for S2, and that it is traveling along the vector connecting 
its current position with the cross indicating the center 
of the nearby system of galaxies, then its inclination from 
the plane of the sky is ~ 75°. Deprojecting according to 
this angle leads to a radius of curvature four times smaller 
and a corresponding IGM density four times larger. 

We can compare the measured IGM densities to those 
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seen in clusters and X-ray bright groups of galaxies where 
the radial density profile of the hot gas in the IGM is 
traced using X-ray data. Gas densities at the center of 
clusters are on the order of 10~ 25 g cm~ 3 (0.1 cm -3 ) 
to 10~ 28 g cm" 3 (10 ~ 4 cm~ 3 ) at a radius of 1 Mpc 
(jVikhlinin et al.l l200l . The IGM density we measure 
near SI is comparable to the density seen in X-ray ob- 
servations at a similar radius in the NGC 5044, NGC 533 , 
and ESO 5520200 groups from lGastaldello et all (|2007l ). 
Those groups have higher velocity dispersions and may 
be dynamically older than the group containing SI. 

Numerical simulations of the evolution of the IGM 
predict WHIM gas in the local universe with a density 
distribution that is broadly peaked at overdensities of 
10 — 20 and densities in filaments of 10 — 100 relative 
to the critical density (1.0 x 10~ 29 g cm~ 3 for H n = 
73 km s" 1 Mpc -1 ) (|Cen fc Ostrikerl [20061 : iDolae et all 
2006). The gas that S2 is tracing has an overdensity of 
~ 90. Considering that this source is 2 Mpc from the 
center of a large group of galaxies it could be tracing gas 
in a filament or gas that is only tenuously associated with 
the nearby group. 



Radio galaxies with bent jets in groups of galaxies pro- 
vide evidence that intragroup gas exists in these systems 
even when the IGM is not X-ray bright. The densities 
derived in this paper indicate that groups of galaxies and 
their surroundings are likely to contain significant reser- 
voirs of baryons. Further work is underway with a larger 
number of sources and will attempt to quantify the con- 
tribution of the intragroup medium to the baryon density 
in the local universe. 
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TABLE 1 

Source Information 



Source 


a 

(J2000) 


S 

(J2000) 


z 


Z/1440 
W Hz" 1 


h 

(arcsec) 


R 

(arcsec) 


Vgal 

( km s 1 ) 


p 

L mm, jet 

(dynes cm~ 2 ) 


llGM 

( cm -3 ) 


SI 
S2 


12 49 42.2 
08 53 31.9 


30 38 38 
23 24 00 


0.194 
0.306 


1.61 x 10 25 
1.88 x 10 25 


4.1 

5.8 


13±1 
23±2 


250+ 20 
zou -100 

570 ± 60 


2 X 10" 11 
7 x 10" 12 


4 ± ]+! X 10~ 3 
9± 3ts° x 10" 4 



Note. — L1440 is from Blanton ct al. (2001). 
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Right Ascension (J 20 00) 



Fig. 1. — 610 MHz GMRT radio contours overlaid on SDSS r band optical images, SI (top) and S2 (bottom). Lowest contour levels 
are 1 mjy beam -1 (SI) and 1.3 mjy beam" 1 (S2) and increase by y/2. Beam sizes are 4.3" X 3.3" for SI and 5.8" X 4.9" for S2. In the 
panels on the left the radius of curvature is illustrated. In the center panels, crosses indicate the averaged group center and galaxies with 
spectroscopic redshifts are circled twice while photometric redshifts are circled once. To the right, the redshift histograms have a dashed 
line indicating the group redshift and a dotted line at the redshift of the bent-double radio source. 



